Distinctive features in the quantum sub-dynamics of bipartite systems 
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Distinctive features in the sub-dynamics of an interacting bipartite system are explored. The 
sub-dynamic Heisenberg operators are introduced in a novel way leading to the quasi-particle-like 
concept in the respective subsystems. This formulation is illustrated using the known exact solutions 
of the generic (Jaynes - Cummings) model of interacting bipartite system of a spin 1/2 particle and a 
single mode electromagnetic field (more generally, boson field) . The striking effects of interaction and 
entanglement are reflected in the dressing of the non-interacting number and spin representations, 
displaying quasi-particle-like properties. This development may be at the heart of the physics of 
back - action discussed in the recent literature. 
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In condensed matter physics, the low- lying excitations 
near the ground state - single particle and collective exci- 
tations - are described in simple terms using the concept 
of quasi-particles (QP). They describe the renormalized 
energy and wave-function due to interaction, but retain 
the pristine nature of the non- interacting particles (e.g., 
fermion nature of the non-interacting electron persists in 
the QP regime). This concept has been widely used in 
more complex situations than their first recognition in 
describing interacting electrons and phonons in metals. 
Many important physical properties of the system at low 
temperatures could be understood in this way. Recently, 
the QP concept has been invoked in exotic systems such 
as Bose-Einstein condensation in two dimensions, quan- 
tum Hall effect, and strongly correlated electron and spin 
systems. It may also be mentioned that the idea of back- 
flow also became evident in this picture. For a compre- 
hensive review of all these, one may consult [l[. 

With the advent of recent advances in experimental 
techniques an entirely different class of problems involv- 
ing time-dependent phenomena, e.g., interacting bipar- 
tite system of atoms and radiation 0, have opened 
up a new class of issues. Even though the total system 
evolves in time unitarily, the time evolution of the subsys- 
tems are quite different [H • By introducing sub-dynamic 
Heisenberg operators, we find a novel dynamical descrip- 
tion of the subsystem time evolution having certain dis- 
tinctive features, which appear to be similar to QP. In 
view of this observation, for want of a better nomencla- 
ture, we suggest quasi-particle-like (QPL) for describing 
these time evolutions. The purpose of the present paper 
is to study the QPL sub-dynamics of both parts of an 
interacting bipartite system on equal footing. In partic- 
ular, the new Heisenberg operators define QPL, but do 
not retain the pristine nature of the non-interacting par- 
ticles, except at certain intervals of time, unlike the QP. 
As a byproduct, this formulation enables us to examine 
the complimentary back-action effects, like the back-flow 



concept mentioned above in the QP context. But unlike 
in the condensed matter systems, QPL describes time- 
dependent features that occur in certain regions of time. 
These dynamical properties of systems, including back- 
action, are being pursued by various experimental groups 
around the world. 

Quantum-limited solid-state devices have shown that 
the model of an atom in a cavity, bathed in a laser field, 
serves as a generic theoretical model for describing these 
systems. The Jaynes-Cummings model (JCM) [5[ serves 
as an exactly solvable model of such two interacting dis- 
parate systems - a two-level system and an electromag- 
netic field (EMF) (or more generally, a boson field) - 
under certain simplifying assumptions. These are the 
rotating wave approximation and weak coupling of the 
two-level atom with the field. JCM has recently been 
extended to strong coupling regime Q. Inherently, this 
model has the atom entangled with the photon and hence 
the dynamics of the subsystems of both the photon and 
the atom will exhibit consequent features of dynamical 
correlation. Although the sub-dynamics of the photon 
and spin (2-state atom) are generally expected to be 
complex for the interacting system, the Heisenberg sub- 
dynamic operators are shown here to make explicit the 
physical processes, which deviate from the pristine dy- 
namics. This allows the identification of QPL regimes 
resembling the pristine dynamics, with renormalized val- 
ues for photon number and spin. These features emerge 
as novel consequences of the formalism presented here. 
It is worth pointing out the use of JCM in several other 
situations in the recent literature: entangled systems for 
controlling symmetric qubits in trapped ions [7j, linear 
optics methods to generate symmetric qubits |8(, collec- 
tive atomic spin excitations [9|, and stability of atomic 
clocks jlOj . 

We first give a brief description of the theory of sub- 
dynamics in terms of the Kraus operators for a typical 
pair of systems 0). We then construct the Heisenberg 
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representations of typical operators in the two subsys- 
tems. We apply this formalism to JCM and discuss how 
this relates to QPL features. This study leads to further 
implications in more general situations. 

Consider two systems A and B with the Hamiltonian 
(time-independent in this sequel) 



H = Hz 



Hp 



Hai 



(1) 



as a generic interacting composite system. The time 
evolution of the density matrix of this system, p(t), is 
expressed in terms of the unitary evolution operator, 
U{t) = e~ ltH , in the standard way , p(t) = U(t)p{t = 
0)W(t). If Os is a typical operator in the composite 
space in the Schrdinger representation, its Heisenberg 
representation is defined by: Tr[Osp(t)] = Tr[0#(£)p(t = 
0)], H (t) = W(t) O s U(t). Here 'Tr' stands for the trace 
over the whole system. When the initial density matrix is 
of the un-correlated form, p(t = 0) = pa(0) <£> Ps(0), the 
time-evolution of these respective density matrices define 
the Kraus representations, with the following expressions 
for the density matrices for the sub-systems: 

PA (t) =Tr B [p(t)] = ^(*W(o)v/(t), 

V?(t)V 4 (t) = Ia, (2) 
p B (t) =Tr A [p(t)] = Y, W } {t) PA {Ci)W](t) 

3 

E ^/WWiW = !b- (3) 
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Here I a, Ib are unit operators in the respective Hilbert 
spaces of the systems A and B. The operators {Vi(t)}, 
{Wj (t)} are expressed in terms of the matrix elements ( in 
the corresponding subsystem spaces) of the unitary time 
evolution operator U(t) and will not be made explicit 
at this juncture. The second set of conditions in Eqs. 
represent the unit traces of the density matrices. 
Consider a typical Schrodinger operator, as , in the 
A-space, which in the composite space is represented 
by a s ® Ib- Then, Tr[a s <g> I B U(t)p(0)W (t)] = 
TtaTtb [a s ® I B U(t)p(0)W(t)} = Tr s [a s p A (t)} = 
TiB[a(t)pA(0)], and thus the effective time-dependent op- 
erator in the space of A is 



&(*) = E V?®*sVi(t). 



(4) 



Similarly a typical time-dependent operator in the B- 
space has the form 



b(t) = Y,W](t)b s W j (t). 



(5) 



Because these transformations are not unitary, the pris- 
tine algebra of the Schrodinger operators is not preserved 
and thus the dynamical effects of interactions are re- 
flected in these modifications. For example, 



(~a(t)) k d k (t), 
a Sl (t)a S2 (t) ^ a sl a S2 (t), 



We now apply this formalism to JCM 0, H, [B[ and illus- 
trate these features in some detail. 

The JCM Hamiltonian and its exact solution La, [a , [11 j: 
We use the convenient notations given in Ref. [1 11 ] with 
units, where the Planck constant h = 1 . We present here 
an alternate description of JCM in terms of the quan- 
tized electromagnetic field (EMF) represented by three 
mutually orthogonal vectors. The direction of propaga- 
tion of photon (say z-direction), electric field perpendic- 
ular to the propagation direction, say y-direction, and 
the magnetic field, transverse to both of these along the 
x-direction. In terms of the creation and annihilation 
operators of the photon field [I3 |: 



E y = i^(uj/2) (a ~ a f ), H x = (a + a 1 ) (6) 

where u> is the photon frequency (laser field). The Hamil- 
tonian of the JCM is then given by: 



H 



E v + H x uo 
ir.\r = — ® I a + y 1 R ® °* 

2 

H — g (H x <E> o~ x + E y <E> a y ), 

i + 1 w o 
= uj (a 1 a + -) (g) I A + — Ir ® (J z 

+g (a (8 a+ + a> (g) a-). 



(7) 



The first term in Eq. ([7|) is the Hamiltonian of the radi- 
ation with frequency uj. The second term is that of the 
atom with the energy difference between its ground and 
excited state ojq and the last term is the dipole interaction 
between the atom and radiation expressed in the rotating 
wave approximation, which leads to dynamical, quantum 
entangled behavior. The second form of the JCM Hamil- 
tonian is the familiar one and opens interpretation in 
more general terms of boson fields in place of EMF. Here 
we employ the standard Pauli matrices and the photon 
destruction and creation operators in the number repre- 
sentation with the standard equal time commutation rule 
[a, of] — Ir 0, EH- This Hamiltonian has a constant of 
the motion, C = a^ o®Ia + \ Ir®cf z = N®Ia + \ Ir®o~ z , 
that enables us to find the normal modes of the JCM. The 
respective unit operators in the radiation and atomic sys- 



tems are Ir = 



E 

ra=0 



|n)(n|, Ia = | T)(t I + I I- The 



time evolution operator in terms of the exact solutions is 
then given by [11] . 



n=0 



+e- fa,t(n+ ^ [«»(*) KT)(n,T I 
-iw n {t) |n+l,|)(n,T | -iw n {t)\n,1)(n + 1,1 Ml 



etc. 



with 

v n (t) 



the photon-atom 



e 1 A "* sin 2 i 



correlation factors 
w n (t) = 



e l A " * cos 2 
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sin 20 n sin A„ t, and 

,2 



tan 9 n = 

1/2 



where 



with detuning frc- 



quency Aw = lo — luq. Consider the initial density 
matrix to be p(0) = p l R ® p\, where p l A denotes 2x2 
density matrix of atoms and p l R = \a)(a\ represents the 

pure state of the laser, with \a) = e I " I / 2 —= |n), 

n=0 

a = \a \ e 1 ^ representing the amplitude and phase of the 
laser beam. The mean number of photons in the beam is 
given by M = \a\ 2 and the phase average of the density 
matrix gives the mixed Poissonian distribution of the 
photons, p(n) = e n! M " . The time-dependent atomic 
marginal density matrix is 

oo 

PA (t) = ]T w NaP \wl a 

W Na = (N\U\a), Wl a = (N\tf\a). (8) 

Similarly the time-dependent radiation marginal density 
matrix is found to be 

PR® = Yl (Pa)s2 S1 V SS2 p r V s \ s , 

8,81 ,S 2 = T,i 



,\U\s 2 ) 



(9) 



It is important to note that the photon sub-dynamics de- 
pends on the initial atomic density matrix and the atomic 
sub-dynamics depends on the initial photon density ma- 
trix. We now express some typical sub-dynamics opera- 
tors in the Heisenberg representation. 

Photons: The annihilation operator in the photon sub- 
dynamics is 



5(t) 



E (pa) 



■-2*1 ^Si S a VsS2 



= e- itu J2[\n)(n + l\V^+TA n (t) 



(10) 



n=0 



+\n-l)(n + l\C n (t) + \n) (n\ D n (t)\ , 



where A n (t) = p\^ (v*(t)v n+ i(t) + w n (t)w n+1 (t) 

P\l (<(*K-l(t) + W n (t)w n -!(t) yfifc) , 

C„{t) = ip\ } (iu n _i (*)<(*) \/n+T- Wn(*X_i(*)\/ra), 
D n (t) = ip' l n (w n -i(t)v n (t) sfn - w n {t)v n -\{t)\/n + l) . 
This is a quasi-annihilation operator of the effective 
photon. The values A n = 1, C n = D n = correspond 
to the usual pristine photon annihilation operator. The 
factor A n represents the suppression of single photon 
annihilation, while C„, D n represent 2-photon and 
no-photon annihilations respectively. The precise nature 
of the photon-QPL will depend on the choice of the 
parameters and on the initial spin state. The standard 
commutation rule, [a, oJ] ~ I R , does not hold for the 
photon-subsystem operators. 



The number operator, used to construct the number 
representation (Fock space) in the non-interacting sys- 
tem, takes the form 

N(t) = (Pa)s 2 si V s \ s a) V S , S2 ? a\t)a(t) 

S,Si ,5 2 
OO 

N(t) = J2 [\n){n\n+\n){n\(p] ]W l(t) - p* u 

n=0 

- ip l n \n+ l)(n| w n (t)v n (t) (11) 
+ ip\\ \n){n+l\w n {t)v n {t)] . 

These relations exhibit how the JCM interaction modi- 
fies the pristine field operators in important ways, thus 
leading to quasi-number representation in the radiation 
subsytem. 

Spin: Using the standard representation of the spin op- 
erators we obtain the quasi-spin components associated 

oo 

with the atom-subspace, &k{t) = E ^Na a 'k(0)WNa' 

s+(t) = ^ (IT>U|Si + (t) + UXT|s 2 + (t) 

+ITXt|S 3 + (*) + UXt|S+(*)) (12) 

= [*-(*)]*. 



with 



stm = e »;ck-i(*) 

n = 

S+(t) = £ »n((K-l(t) 



a* -Jn 



a \/n + 1 
s/ri 



And, 



with 



71 = 

^ Q Vn + l 



**(*) = ITXT|Sf(*) + U>UI^(t) 
+ IT)ai^W + IIXTI^W, 



(13) 



slit) 
slit) 



1-2 £ v(n)w 2 n (t) 



1-2 E P(™ + !) ™n(*) 



-2« E P(") *(<KW ^-^r = ST(t) 



Here again the time-dependent coefficients display the 
contributions from the JCM interactions and they de- 
pend on the initial photon state. The algebra of these 
operators is different from that of the pristine counter- 
parts as can be seen directly. It is emphasized that these 
transformations are not unitary, the pristine algebra of 
the Schrodinger operators is not preserved and thus the 
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dynamical effects of interactions are reflected in these 
modifications. 

Results: Explanation of QPL features. 

The sub-dynamics of the photon and spin are compli- 
cated due to the presence of new photon and spin pro- 
cesses as seen in Eqs. (fl"0|) . (fl2l and (IT3)) resulting from 
JCM. The photon QPL properties will emerge if its sub- 
dynamics is dominated by A n in Eq. IjlOp and the other 
two terms negligible. Similarly the spin QPL proper- 
ties correspond to a regime, where the first two terms in 
Eq. (fT3|) dominate, while the last two terms are negligible. 
Under these circumstances, the sub-dynamics appears to 
admit the interpretation in terms of fractional numbers of 
photons and fractional spins. Such QPL regimes emerge 
when detuning dominates the Rabi term of the coherence 
factor X n , i.e., Aoj » 2g \Jn + 1. Then the resulting dy- 
namics of QPL photons and spins fall into collapse-revival 
structure of a large detuned JCM. In Fig. 1, we display 
how far the operator structures of the photon, Eqs. (fT0| . 
(jlip . and the spin, Eq. (|13p deviate from their pristine 
non-interacting versions in the subystem for various val- 
ues of the system parameters. In the non-interacting 
case, (a\a(t)\a) — a, (a\N(t)\a) = \a\ 2 and two sets 
of curves in Fig. (la) show the effects of interaction on 
these values as a function of gt. In Fig. (lb) we display 
the eigenvalues of <j z {t) vs. gt. These eigenvalues are 
given by (OS% ± (ds) z , with the off-set from zero is given 
by (OS)* = (Sf(t) + s^*))/ 2 : an d the dispersion is given by 
(DS)z = j^ < s fw- s iM> ) 2 + |gg | 2 j> 1/2 . In the non-interacting 

case, these are +1 and -1, with off-set equal to 0. The 
curves in Fig. (lb) show these and, in addition, display 
the collapse and revivals. The collapses mimic constant 
values different from ±1, depending on the parameters 
and the revivals occur at different times depending on 
the detuning. Also, at half the collapse time, the pho- 
ton figure shows a minimum, and at peaks of revivals, 
the photons exhibit a maximum. We must point out 
that similar results are obtained when other choices of 
parameters were made. The off-set being small and the 
dispersion different from +1,-1 and fractional, are the 
results of the non-unitary evolution governing the sub- 
dynamics. A more detailed investigation of these features 
will be explored in another communication. 

The back action phenomena may also be interpreted 
in another complementary way. We have already shown 
how the pristine representations of spin and radiation 
get modified due to interaction and entanglement. We 
now exhibit this in terms of the constant of motion, C. 
Upon taking the expectation value of this operator over 
the total density matrix of the system, we obtain, (n) + 

3 (Ptt ~ PlO = + 3 foW) ' where ^ and 

are given explicitly by Eqs. (fTTj) and (fl3j) respectively. 

The back action on the photon number (N(t)) is such 

that it compensates the collapse and revivals that occur 

in (a z (t)). Thus the photons and the atom act in tandem 



in such a way as to maintain the constant of the motion. 
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FIG. 1: (a) Magnitude of (a\a(t)\a) and (a\N(t)\a) vs. gt. The set of 
JCM parameters used are, Atom: p^(0) — 1, all others zero; Photon: 
jV mean = | Q | 2 = 10; Interaction: Red: Aw/g = 7.5, Blue: Aw/g = 10 
and Green: Au) / g — 20. Fig. 1 (b) displays the offsets and eigenvalues 
of a z {t) vs gt for 7V moa „ = 10. 



In conclusion, we have displayed here the complemen- 
tary dressed behaviors of the sub-dynamics of photons 
(equivalently, the EMF) and atom in the JCM model of 
interacting fields. We believe that this type of behavior 
holds generally in the non-unitary sub-dynamic evolution 
of any bipartite system. The details will differ depending 
on the type of interaction between the two systems. 
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